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ABSTRACT The HiRes algorithm has been used to analyze the far- 
infrared emission of shocked gas and dust in supernova remnants. In the 
case of supernova remnant IC 443, we find a very good match between the 
resolved features in the deconvolved images and the emissions of shocked 
gas mapped in other wavelengths (hues of H 2 , CO, HCO+, and HI). 
Dust emission is also found to be surrounding hot bubbles of supernova 
remnants which are seen in soft X-ray maps. Optical spectroscopy on 
the emission of the shocked gas suggests a close correlation between the 
FIR color and local shock speed, which is a strong function of the ambi- 
ent (preshock) gas density. These provide a potentially effective way to 
identify regions of strong shock interaction, and thus facilitate studies of 
kinematics and energetics in the interstellar medium. 


INTRODUCTION 

Supernova remnants (SNRs) can be bright in the mid- and far-infrared ei- 
ther because of noil-thermal radiative processes internally, or interactions with 
the interstellar medium (ISM) surrounding them. Several recent studies have 
surveyed a large number of Galactic SNRs in the IRAS All Sky Survey, and iden- 
tified a substantial fraction of them as FIR sources despite considerable confusion 
problems due to the Galactic plane (Arendt 1989; Saken et. al. 1992). General 
correlations between IR and radio, and/or between IR and X-ray emissions are 
cited, even though the resolution that IRAS data provide is often insufficient 
to study the IR morphology. In addition, it is found on a statistical basis that 
older remnants are both colder and more massive. These results suggest that 
the emission measured by IRAS come mostly from shock interactions with dense 
gas and dust. However, the actual mechanism of interaction, and its dependence 
on local conditions largely remain to be explored. 

An alternative approach to the statistical analysis is trying to better resolve 
the IR emission around a few extended SNRs, and compare the images in the 
four IRAS bands with the physical conditions in the surrounding medium studied 
in other wavelengths. In this way, one can hopefully address the origin of the 
infrared emission as well as the energy balance in the evolution of the remnants. 
The HiRes algorithm (Aumann et. al. 1990) becomes particularly useful in this 
regard, because it provides not only a much better linear resolution, but also 
a convenient way to examine infrared color variation through beam matching. 
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In this contribution, I present the preliminary results of a HiRes analysis to 
illustrate some of the possibilities that this approach can provide. 

DATA ANALYSIS 

The FIR emission of supernova remnant IC 443 is detected in all four IRAS 
bands and its ring-like structure has been well identified (Huang et. al. 1986; 
Braun & Strom 1986; Arendt 1989). However, since the angular scale of the 
entire region is about 20 ; , the ~ 4' resolution of the coadded images is inadequate 
to analyze the infrared morphology. As a result, only the integrated fluxes have 
been measured in the previous studies. 

Since both the absolute flux level and the signal-to-nise ratio in the 60 and 
100 pm images of SNRs are higher than that in 12 and 25 pm, it is justified to 
perform more HiRes iterations on the longer wavelength images. In practice, I 
used the survey data as input, and adopted the HiRes 12, 25, 60 and 100 pm 
images at 5, 5, 10, and 20 iterations, respectively. Although these do not provide 
the best possible resolution at all four bands, they do give more comparable 
beam sizes. More importantly, they are physically meaningful because of the 
different signal-to-noise ratios in the four bands. The region selected to perform 
the deconvolution is a 2° field centered on IC 443. 

The variation of the beam from the field center to the sides is rather small in 
this case, so I adopted measurements of simulated beam at the image center to 
be the representative one. The next step was to examine the resolution maps at 
each band, and try to match the resolution based on the “least common denomi- 
nator ’. Four different two-dimensional smoothing kernels were used to convolve 
the maps into nearly identical beam sizes and shapes, while preserving the best 
possible resolution. The parameters of the kernel used, and the measurements of 
input and output beams are listed in Table 1. These measurements were based 
on 2-D least-square Gaussian fits. 


Table 1. Beam measurements" 


Before Convolution 6 Gaussian Kernel" 


After Convolution 6 


Band(0 d 

a(maj) 

a(min) 

sa 

ar 

a(maj) 

a(min) 

b 1 (5) 

6.58 

2.11 

1.65 

0.920 

7.55 

4.45 

b2(5) 

7.16 

2.22 

1.62 

0.573 

7.56 

4.45 

b3(10) 

5.34 

3.12 

2.25 

0.589 

7.54 

4.45 

bl(20) 

5.33 

3.97 

2.25 

0.367 

7.55 

4.46 


“ all measurements are in unit of pixel, 1 pixel = 15"; 

6 a(maj), a(min) are half width along the major and minor axes, respectively; 
c sa, ar are the semi-major axis and axis ratio of the 2-D Gaussian; 
d bn(i): where n represents band id, and i is the number of iterations used. 
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From Table 1, we see that the beams in the four IRAS bands are well 
matched after convolution with the selected kernel, with a relatively minor 
penalty of a ~ 30% increase in beamwidth. The beam matching of the ac- 
tual images was confirmed by comparing a few stellar sources present in the 
same field. Figure 1 compares the beam-matched 60 and 100 f-trn images with 
the unprocessed ones. The resolution of the HiRes images is about 1.9' x 1.1', 
a significant improvement over the ~ 4' X 3' beam in pre-HiRes images. The 
biggest advantage of having beam-matched maps in all four bands, however, is 
to make infrared color maps that would be physically meaningful. One can then 
easily compare colors of different regions, provided that the sizes of regions over 
which the flux integrations are made is larger than the beam size. 

COMPARISON WITH OTHER OBSERVATIONS 

Before trying to interpret the physical meaning of the HiRes images and 
color maps, it is necessary to assess the reliability of the enhanced features seen 
in the HiRes deconvolution. In IC 443, this can be best shown in comparisons 
with observations of the same region made in other wavelengths. 

Using millimeter- wave radio telescopes, Dickman et, al. (1992, see also 
Huang et. al. 1986) mapped the broad CO and HCO+ line emission in the sur- 
rounding areas of IC 443. These molecular lines, typically as broad as 60 km s -1 , 
are often identified as the shock interface with nearby molecular clouds. Over- 
lapping the CO and HCO+ maps with the deconvolved IRAS images, we find 
that all broad molecular line regions also have strong FIR emission. These re- 
gions are located mainly along the cj-shaped ring structure in the middle and 
southern parts of the SNR. Some of the previously unresolved FIR features, 
such as the one corresponding to “position F” (Dickman et. al. 1992) near the 
geometric center of IC 443, are clearly isolated FIR emission peaks, similar to 
their counterparts in the molecular hue survey. 

Burton et. al. (1988, 1993) made extensive maps of the (1-0) S(l) line 
emission of molecular hydrogen covering the entire region of IC 443. The near- 
infrared line represents shocked molecular gas at a warm (~ 2000 K) temper- 
ature. Although the H 2 map is at a much higher resolution, the overall cor- 
respondence with the deconvolved IRAS images is striking. This suggests that 
the mid- and far-infrared emission detected by IRAS is mainly from heated dust 
associated with the shocked molecular gas traced by the H 2 line emission. This 
comparison serves two purposes: on the one hand, some emission regions with 
structures poorly resolved in the HiRes images can be studied by referring to 
the NIR image; conversely, the FIR color information extracted from comparing 
different IRAS bands can be used to help interpret the significant variations of 
H 2 emission around this SNR (see below). 

Wang et. al. (1994) mapped the 21 cm neutral hydrogen line emission of 
IC 443 with relatively high resolution (15" as compared to the 50" resolution 
in previous studies). Again, the regions with strong infrared and molecular line 
emission also show up prominently with broad-line H 1. Detailed comparisons 
indicate that the detected H I is mainly the diffuse gas resulting from shock 
dissociation of molecules as the SNR interacts with the surrounding dense clouds. 
However, the kinematic information from HI is incomplete due to the foreground 
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Before Hi Res: 


After Hi Res: 






FIGURE 1 Greyscale maps of the IRAS 60 (upper) and 100 //.m (lower) im- 
ages of IC 443 before and after Hi Res processing. This is a one-degree held 
centered on the supernova remnant; north is at top, and east is to the left. Grey 
levels are linear. Some of the vertical stripes seen in the processed images are 
due to artifacts of the deconvolution algorithm. The HiRes images are “beam 
matched”, such that the effective resolutions are very similar (see text). 
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absorption within the Galactic plane. 

Finally, the IPC aboard Einstein and, more recently, the PSPC aboard 
ROSAT both made large-scale diffuse X-ray emission maps of IC 443 in the 
0.1-3 KeV energy range (Pctre et. al. 1988; Rho & Petre 1993). Comparison 
between the infrared and soft X-ray maps is intriguing because the higher en- 
ergy photons are emitted by the hot interior of the SNR bubble (e.g., Arendt 
et. al. 1991). Indeed, their positional coincidence with the rim of X-ray emit- 
ting regions provides independent confirmation that much of the faint structures 
seen in the HiRes images, especially in the NW region of IC 443, are actual SNR 
shock fronts. 

In summary, the HiRes images with much improved resolution clearly demon- 
strate the relationship between the mid- and far-infrared emission and the SNR 
shock interactions in the ISM. 

INFRARED COLOR VARIATION 

As pointed out above, one of the great advantages for having beam-matched 
IRAS images is to construct color maps of infrared emission, which are important 
to the studies of shock interactions. 

In the case of 1C 443, we find that the FIR color changes drastically from 
one part of the SNR shell to the other, presumably tracing the density variation 
in the ambient medium. In the NE part of the SNR where the optical filaments 
are prominent, there is a low S(l'2)/S(25) and a high S(60)/S(100) flux density 
ratio as compared to the rest of the FIR bright regions (SE and SW), where 
molecular line emissions are the strongest. 

Theoretically, the energy spectrum of the dust emission in interstellar shocks 
depends critically on the average preshock density, n 0 . Where n 0 is low, the 
shock speed is high and so is the heating rate. The corresponding FIR (60 and 
100 fan ) color then reflects a higher dust temperature. Conversely, lower dust 
temperatures are found where the preshock density is high. In IC 443, the latter 
case is typically regions dominated by molecular gas. 

We have also obtained spectroscopic (echelle) data of the optical filaments 
at a velocity resolution of about 8 km s _1 . Since the curvature of the shock 
front always tends to broaden the observed linewidtli, it can be shown that 
the FWHP (full-width at zero power) of the emission lines of any particular 
region is a direct indicator of the shock speed at that position. And it is so 
regardless of the actual angle between the shock velocity vector and the line of 
sight (Dopita 1979). Adopting this assumption, we find in the NE where the 
S(60)/S(100) color ratio is high, the shock speed is about twice that of the SE 
and SW regions (120 km s -1 versus 60 km s -1 ). The optical data thus confirm 
the idea that the SNR has encountered very dense gas in the southern part, 
resulting in faster cooling and a dramatic slow-down in the shock propagation. 

The detailed morphology of the regions with different FIR colors can also be 
quite different. Referring to the high resolution H 2 line images, we find that the 
NE region shows thin, multi-layered filamentary structures; while the SE and 
SW regions show more pronounced, ripple-like features (Burton et. al. 1993). 
These difference could be due to the changing characteristics of the shock as its 
speed drops substantially (e.g., Wang & Scoville 1992). 
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DISCUSSION 

One other thing that is clear from the color maps is that an enhanced ratio 
of S(60)/S(I00) is usually accompanied by a decreased S(12)/S(25). This is 
also true for integrated fluxes of a large sample of SNRs (e.g., Arendt 1989). It 
reminds one of a similar trend among galaxies (Helou 1986), which is accounted 
for in terms of changing conditions in the ionizing radiation (Desert 1986). In 
the case of SNR shocks, however, the explanation maybe somewhat different. It 
is likely that the color ratio change among different shocked regions in a SNR 
is directly related to the amount of very small grains (VSG) that survive at the 
shock front. Faster shocks in the diffuse gas (with warmer 60-100 firn emission) 
tend to destroy most of the VSGs which are responsible for the emission at the 
shortest (12 fim) wavelength, resulting in a lower S(12)/S(25) ratio. 

For interstellar shocks in dense gas, line emissions could potentially alter the 
colors measured using IRAS flux density ratios. Burton et. al. (1990) detected 
a strong [01] A 63 pm line at one of the molecular emission peaks in 1C 443, 
and estimated that it may contribute substantially to the local broad-band 60 
fim flux density. If this is generally the case, then the IRAS data, even with 
the beam-matched images, cannot be used to derive effective dust temperatures 
directly. More FIR spectroscopic observations are therefore needed to help better 
address the real properties of dust emission in the shocked ISM. 
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